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This proceeding summarizes the Higgs recoil mass and Higgs-strahlung cross-section 
study done for the Letter of Intent (LOI) of the International Large Detector (ILD) 
Concept. Assuming Mh = 120 GeV, working at y/a — 250 GeV with beam parameters 
RDR250 and beam polarization P{e~ , e^) = (—80%, +30%), this full simulation study 
predicts that, the ILD detector can achieve 37 MeV precision on the Mh measurement 
^1^' and 3.3% on the cross-section measurement from the ZH — >■ fx^ yT X channel, while 

^ ! 83 MeV and 4.9% from the ZH -^ e+e"X channel, if we have 250 fb"^ integrated 

luminosity. 
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>■ ' 1 Introduction 

O/^ ' Experimental conditions at the proposed International Linear Collider [T] (ILC) provide 
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an ideal environment for the precision studies of Higgs production own to the unparal- 
leled cleanliness and vifcU-defined initial conditions. We can use the Higgs-strahlung process 
(e^e^ — > Z* — > ZH), which is a major Higgs production mechanism at the ILC, to precisely 



^— s ' measure the Higgs mass {Mh) independent of its decay modes, using the mass recoiling to 

the Z boson, with Z — )■ ^i^ pL~ or e+e". The major advantage of this recoil mass method is 
we do not need any information about how the Higgs decays. To calculate the Mh, we only 
need the Z boson reconstructed from the lepton pair it decays with the precisely known cen- 

^^ ' ter of mass energy {\/s). We can also precisely determine the Higgs-strahlung cross-section 

H \ (cr) and therefore the coupling strength at the HZZ vertex, without any bias from the Higgs 

decay assumptions. 

This proceeding summarizes the Higgs recoil mass and Higgs-strahlung cross-section 
study done for the Letter of Intent [4] (LOI) of the International Large Detector [2] (ILD) 
Concept. It is a most realistic full simulation study up to now based on the ILC design issued 
in the ILC Reference Design Report [3] (RDR) and the ILD detector concept. Thus, it has 
contemplated all the major effects we have thought of from both the accelerator and the 
detector. It assumes Mh = 120 GeV, the ILC works at ^/s = 250 GeV, and an integrated 
luminosity of 250 fb~^ (equivalent to four years of data taking). The study also assumes 
two options of beam longitudinal polarization, e£ej.: (e~ : —80%, e"'" : -1-30%) and e^e^^: 
(e~ : -1-80%, e^ : —30%), to quantitatively estimate the benefits from the beam polarization. 
The results of this study are published in the ILD LOI, with an accompany physics 
note [5] summarizing the methods. For the complete documentation of this study, please 
refer to my Ph.D thesis [6]. 
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2 Production 

The event generation is centrally done at SLAC [7] using WHIZARD [8] vl.40, based on 
a detailed beam simulation using GUINEA-PIG [H] with beam parameters given in ILG 
RDR, in which the beam energy spread are 0.28% and 0.18% for electron and positron 
beams, respectively. Both the Initial-State Radiation (ISR) and Final-State Radiation (FSR) 
are included in the event generation. The full simulation of the ILD detector and the 
reconstruction are done at DESY and KEK using software package ILCSoft [THj vOl-06. 

Tabic [T] give the cross-sections of signal and major backgrounds considered in this study. 
We categorize the backgrounds according to the final states. For instance, the ee background 
consists of e'^e^ — >■ Zj^* — >■ e^e~, the fivfiv background is mostly the WW — )■ fivfiv, but 
also ZZ — > i^i^iyv, while the ee// background comes from all the possible intermediate states 
e+e- -^ Z/"f*Z/"t*. 
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Cross- 


Section 




eZ^H 


ei?.4 


fifiX 


11.67 fb 


7.87 fb 


fifi 


10.44 pb 


8.12 pb 


TT 


6213.22 fb 


4850.05 fb 


fiiyfiv 


481.68 fb 


52.37 fb 


^^^^ff 


1196.79 fb 


1130.01 fb 



Process 


Cross- 


Section 




eZCfl 


^R^L 


eeX 


12.55 fb 


8.43 fb 


ee 


17.30 nb 


17.30 nb 


TT 


6213.22 fb 


4814.46 fb 


evev 


648.51 fb 


107.88 fb 


ee// 


4250.58 fb 


4135.97 rb 



Table 1: Gross-sections with the two options of beam polarization e^e 
signal is indicated by bold face letters. 
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3 Event Selection 

The first step in the event selection is the identification of leptonically decaying Z bosons. 
We require the lepton tracks to be well measured by removing tracks with large uncertainties 
on the reconstructed momentum. The efficiency to identify a pair of Icptons from the decay 
of a Z boson is 95.4 % for Z -^ n+fj,' and 98.8 % for Z -^ e+e". 

The second step in the event selection is to suppress the backgrounds who also have a 
pair of leptons identified. We can suppress most of them using some common kinematic 
requirements. These requirements are, the invariant mass of the lepton pair (Mdi) is around 
the Z boson mass, a large transverse momentum (Ptcu) of the lepton pair, a not-so- forward 
polar angle (cos 6di) distribution of the lepton pair system, the pair of leptons is not back- 
to-back (acoplanarity) , and a restrict on the recoil mass (Mrecou) window around the Higgs 
mass. 

There are irreducible backgrounds like lepton pairs (^/i or ee), WW and ZZ after the 
requirements above. 

For the irreducible lepton pair background, we have found another method named Pt 
balance. The idea is to identify the energetic ISR photons, and ask them to balance the Pt of 
the lepton pair. We can observe a great Pt balance from the lepton pair background, while 
not from the ZH signal. This is because, the irreducibility of the lepton pair backgrounds 
is owen to the fact that the ISR energy loss reduces the effective -^/s, thus migrates them to 
the acceptable window of our common kinematic requirements. Together with a veto of the 



LCWS/ILC2010 



photon conversions by requiring a large polar angle difference between the pair of leptons, 
we can totally remove the lepton pair backgrounds. 

For the irreducible WW and ZZ backgrounds, we employ a Likelihood method for 
the further suppression. We use variables Mdi, Prdi, cos 6 di and acoUinearity to build the 
templates for Likelihood calculation, and optimize the cut on the Likelihood according to 
the maximum of the significance. The Likelihood further suppression rejects the WW and 
ZZ backgrounds further by one half. 

With this event selection scenario, we can keep an efficiency of around 55 % with a signal 
over background ratio of about 4 under the recoil mass peak for the ZH — > fi/^iX channel. 

4 Fitting Methods 

After the event selection, the remaining spectrum consists of the ZH signal (S) and the 
irreducible backgrounds (B). Now we need to build a composed model Fm{x; Mff,Ns) = 
Ns ■ Fs{x; Mh) + Nb -FBix), and fit it to the remaining spectrum to extract our results. 
In the composed model, the Mjj and Ns are the free parameters we need to determine. M^ 
is the Higgs mass, and the Ns is for the calculation of the cross-section. 

We choose the polynomial function to describe the shape of the background. For the 
signal, I have studied three functions in my thesis. They arc: 

• Use Gaussian to describe the Peak and adding an Exponential to complement the Tail 
(GPET). It is function modified from previous contributions [TH[TT]. My modification 
makes both itself and its first derivative being continuous. 

• An implementation of the Kernel Estimation |13j using Gaussian kernel. 

• A Physics Motivated function first developed in my thesis. The development started 
from Yokoya-Chen's Bcamstrahlung approximation [TH (TSl HH HZ] , analytically con- 
voluted with the ISR approximation [18j . After a numerical convolution with the 
Gaussian function, I propagated it to the recoil mass. With this Physics Motivated 
function, and known the beam parameters, I can essentially predict the generator level 
ZH spectrum. 
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Figure 1: Fits using the Physics Motivated function to describe the signal, for ZH 
HfiX (Left) and eeX (Right) channels, taking beam polarization e£e^ as examples. 



LCWS/ILC2010 



Figure [T] shows the fits using the Physics Motivated function to describe the signal, for 
both the ZH — > fjifiX and eeX channels, taking beam polarization e£e^ as examples. The 
precisions derived from the fits are 37 MeV on the Mh measurement and 3.3 % on the a 
measurement from the fjifiX channel, while 83 MeV and 4.9 % from the eeX channel. 



5 Bremsstrahlung Recovery 

The results derived above shows the eeX channel gives almost twice the errors on both 
the Mh and the a measurements, than those from the /i/iX channel. This is because 
the electrons suffer heavily from the Bremsstrahlung energy loss, reminding that the ILD 
detector has about 4 % Xq material budget before TPC. 
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Figure 2: Left: Comparison of the recoil mass distributions with and without 
Bremsstrahlung recovery. Right: Fit of the eeX channel with Bremsstrahlung recovery 
corresponding to that in Figure [TJ 

A dedicated algorithm |19| is employed in this study to identify the Bremsstrahlung 
photons and merge them into the Z boson. Figure [2] (Left) shows a comparison of the 
Higgs recoil mass spectra with and without Bremsstrahlung recovery. We can observe the 
Bremsstrahlung recovery significantly increases the statistics while degrades the mass reso- 
lution because of the worse ECAL energy resolution for low energy photon reconstruction. 

Figure [2] (Right) shows the fit of the eeX channel with Bremsstrahlung recovery cor- 
responding to that in Figure [T] The results with Bremsstrahlung recovery are 73 MeV on 
the Mh measurement (~ 10 % improvement), and 3.9 % on the a measurement. With the 
Bremsstrahlung recovery, the accuracy of the Mh measurement is still worse than that of 
the ^b^X channel by a factor of two. However, the accuracy of the cross section measurement 
becomes similar to that of the fi^X channel, since it is less sensitive to the mass resolution. 



6 Results 

Table [5] summarizes the resulting precisions of the measurements for both beam polarization 
^~L^R and Cj^ej^; where those of the eeX have Bremsstrahlung recovered, and the "merged" 
means the results by combining the two Icptonic channels. 
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Pol. 


Ch. 


Mh (GcV) 


a(fb) 


eZ4 


Ai/iX 


120.001±0.037 


11.63±0.39(3.35%) 


eeX 


119.997±0.073 


12.52±0.49(3.91%) 


merged 


120.006±0.033 


12.02±0.31(2.54%) 


^R^t 


^fiX 


119.997±0.040 


7.82 ±0.28(3.58%) 


eeX 


120.003±0.081 


8.41 ±0.36(4.28%) 


merged 


120.005±0.035 


8.09 ±0.22(2.73%) 



Table 2: Resulting Higgs mass Mh and cross section a, of the ^jiX channel and eeX 
channel with Bremsstrahlung recovery. The "merged" means the results by combining the 
two leptonic channels. 



7 Conclusion 



The results in Table [5] shows the best results of a single leptonic channel come from ZH -^ 
jifiX. Even including the Bremsstrahlung recovery, the eeX channel gives worse results 
than the ^^X channel by a factor of two in the Mh measurement, and by a factor of 1.5 
in the a measurement. This is because of the worse energy resolution of the soft photon 
measurement by the ECAL, and larger backgrounds related to eeX channel. 

The results with beam polarization e^e^ are better than those of the e]^e^ by about 
10 %. The reasons are that, although the polarization e]^ej suppresses the WW background, 
the cross sections of the Higgs-strahlung process are smaller by about 20 % compared to that 
of the polarization e£e^. At the same time, the methods we have developed are efficient 
enough for the suppression of the WW background. Hence the polarization e^e^ gives 
better results. 

We may also conclude that there are two accelerator effects, the beam energy spread 
and the Beamstrahlung, have great impact on Higgs recoil mass measurement. The beam 
energy spread increases the width of the recoil mass peak, thus reduces the accuracy of the 
mass measurement. The Beamstrahlung largely reduces the effective statistics on the recoil 
mass peak. 
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Figure 3: The Higgs recoil mass distributions in the /x/iX channel (Left) and eeX channel 
(Right), comparison of those before and after detector simulation. 
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Figure [3] compares the Higgs recoil mass distributions before (generator level) and after 
(reconstructed) full detector simulation and reconstruction for fXfiX channel (Left) and eeX 
channel (Right). The distribution in the generator level has the accelerator effects imposed, 
while that after reconstruction has the detector effects added. 

For the recoil mass distributions of /i/iX-channel, a fit to the left side of the maximum 
using a Gaussian function gives the mass resolutions to be 560 MeV in the generator level, 
and 650 MeV after full detector simulation. The detector response leads to a broadening of 
the recoil mass maximum from 560 MeV to 650 MeV. The contribution from the uncertainty 
of detector response is therefore estimated to be 330 MeV. This observation indicates that 
the dominant contribution to the observed width of the ^fiX recoil mass distribution arises 
from the incoming beams rather than the response of the ILD detector. 

8 Outlook 

All what we have intensively studied in the past is how to reduce the statistical errors on 
these two measurements. A serious study of the systematic errors is apparently missing in 
the context. 

In the Higgs recoil mass measurement, the systematic biases appear as the difference 
between the measured value and the true value. Observables correlated with the Higgs 
recoil mass are the center of mass energy and the momenta of the pair of leptons. Imper- 
fect knowledge of these variables could introduce systematic bias to the Higgs recoil mass 
measurement. To control them, a well studied reference reaction is needed. 
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Figure 4: Fit of the Invariant Mass (Left) and Recoil Mass (Right) of the Z 
the ZZ process using Kernel Estimation, with polarization mode e 
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13 MeV is obtained for the invariant mass measurement, 
measurement. 



An accuracy of 
and 28 MeV for the recoil mass 



L'^B- 



The ZZ process is an excellent choice, with the Mz precisely known to a precision of 
2 MeV [5D] . The ZZ process has a similar scenario as the Higgs-strahlung process with one 
Z boson decays to a pair of muons or electrons, and the recoil mass of the Higgs replaced 
by that of the other Z boson. 

By measurement of the invariant mass of the Z — > l^l^, the tracking momentum 
measurement can be calibrated. At ^/s = 250 GeV, assuming Mh = 120 GeV, the 
ZZ -^ iniX/eeX has about 40 times larger cross section than that of the ZH -^ fif^iX/eeX. 
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With this much larger statistics the Z mass can be measured to a precision of 13 MeV, using 
channel ZZ — )- n^X, see Figure H] (Left) for the fit. 

The Z recoil mass of the ZZ process can be used to determine and control the center of 
mass energy and the radiative effects. The Z recoil mass could be determined to a statistical 
precision of 28 MeV, using channel ZZ — ?► fifiX, see the fit in Figure |4] (Right). With this 
small statistical error, the knowledge of the center of mass energy and the radiative effects 
could be validated precisely. 

On the cross-section measurement, the major systematic error comes from the uncer- 
tainty of the efficiency. The procedure to measure the uncertainty of the efficiency is to vary 
the physics assumptions together with the various background rejection methods to estimate 
the dependences and covariances between them. Reminding we applied a very complicated 
background rejection method. The complicated background rejection mandatorily intro- 
duces more difficulties in the measurement of the uncertainty of the efficiency. 

For the cross section measurement, we found that the statistical error does not request 
such a high suppression of the background as for the mass measurement. By removing 
the Likelihood method and remaining only cuts on some basic variables Mdi, Prdi, and Pt 
balance, the statistical error on the cross section measurement only increases by about 10 % 
on average. This means a similar statistical error could be obtained with much less sources 
of systematic error due to the background rejection methods. 
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